US  Army  Corps 
of  Engineers® 

Engineer  Research  and 
Development  Center 


INNOVATIVE  SOLUTIONS 
for  a  safer,  better  world 


A  Creep  Model  for  High-Density  Snow 

Robert  B.  Haehnel  April  2017 


Approved  for  public  release;  distribution  is  unlimited. 


The  U.S.  Army  Engineer  Research  and  Development  Center  (ERDC)  solves 
the  nation's  toughest  engineering  and  environmental  challenges.  ERDC  develops 
innovative  solutions  in  civil  and  military  engineering,  geospatial  sciences,  water 
resources,  and  environmental  sciences  for  the  Army,  the  Department  of  Defense, 
civilian  agencies,  and  our  nation's  public  good.  Find  out  more  at  www.erdc.usace.armv.mil. 


To  search  for  other  technical  reports  published  by  ERDC,  visit  the  ERDC  online  library 
at  http://acwc.sdp.sirsi.net/client/default. 


ERDC/CRREL  TR-17-7 
April  2017 


A  Creep  Model  for  High-Density  Snow 


Robert  B.  Haehnel 

U.S.  Army  Engineer  Research  and  Development  Center  (ERDC) 
Cold  Regions  Research  and  Engineering  Laboratory  (CRREL) 
12  Lyme  Road 
Hanover,  NH  03755-1290 


Final  Report 


Approved  for  public  release;  distribution  is  unlimited. 


Prepared  for  Smithsonian  Astrophysical  Observatory 
60  Garden  Street 
Cambridge,  MA  02138 

under  Project  SA0-GLT-16-01,  “Greenland  Telescope  Snow  Foundation  Performance 
and  Construction” 


ERDC/CRREL  TR-17-7 


ii 


Abstract 

This  effort  modifies  Mellor  and  Smith's  creep  model  for  dense  snow  to 
conform  to  the  more  general  creep  power  law  form  (Glen's  creep  law  for 
ice  is  a  special  case  of  this  power  law) .  The  present  study  used  this  general 
form  as  the  basis  for  developing  two  creep  models:  one  to  describe  the  pri¬ 
mary  creep  and  another  that  models  secondary  creep  response  as  a  func¬ 
tion  of  applied  stress,  snow  density,  and  temperature.  Available  laboratory 
data  helped  to  determine  the  parameters  for  these  models.  These  models 
were  recast  into  a  form  compatible  with  the  ABAQUS  finite  element  soft¬ 
ware  to  allow  simulation  of  settlement  of  a  snow  foundation  subj  ected  to  a 
long-term  load. 

The  secondary  creep  model  was  validated  against  the  observed  settlement 
of  a  spread  foundation  on  snow  at  Camp  Century,  Greenland.  Though 
there  were  not  suitable  field  data  to  validate  the  primary  creep  model, 
comparison  of  that  model  to  results  obtained  for  the  secondary  creep 
model  demonstrate  that  the  primary  creep  model  behaved  as  expected. 
Additional  work  is  necessary  to  merge  the  primary  and  secondary  creep 
models  and  to  further  validate  the  primary  creep  response  of  the  model 
presented  here. 


DISCLAIMER:  The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication,  or  promotional  purposes.  Ci¬ 
tation  of  trade  names  does  not  constitute  an  official  endorsement  or  approval  of  the  use  of  such  commercial  products. 
All  product  names  and  trademarks  cited  are  the  property  of  their  respective  owners.  The  findings  of  this  report  are  not  to 
be  construed  as  an  official  Department  of  the  Army  position  unless  so  designated  by  other  authorized  documents. 

DESTROY  THIS  REPORT  WHEN  NO  LONGER  NEEDED.  DO  NOT  RETURN  IT  TO  THE  ORIGINATOR. 
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1  Introduction 

1.1  Background 

The  Greenland  Telescope  (GLT)  is  a  radio  telescope  that  the  Smithsonian 
Astrophysical  Observatory  proposes  to  install  near  Summit  Station,  Green¬ 
land,  and  that  will  be  founded  on  a  compacted  snow  surface.  The  defor¬ 
mation  of  snow  under  a  constant  load  (creep  deformation,  or  simply 
"creep")  is  dependent  on  both  the  density  and  temperature  of  the  snow. 
Mellor  and  Smith  (1966),  Mellor  (1969),  Reed  (1966, 1974),  Kovacs  (1967), 
Meussen  et  al.  (1999)  and  Chandel  et  al.  (2007)  have  previously  studied 
this  creep  behavior.  Creep  is  normally  characterized  as  proceeding  in  three 
phases  as  depicted  in  Figure  1.  During  the  initial  phase,  termed  primary 
creep,  the  strain  rate,  s,  varies  with  time,  t,  following  roughly  the  trend 
e  oc  t~p  where  p  is  a  creep  exponent.  Following  this  initial  period,  second¬ 
ary  creep  sets  in,  during  which  the  strain  rate  is  constant  with  time.  In  ter¬ 
tiary  creep,  the  strain  rate  once  again  is  time  dependent  following  the 
trend  of  s  oc  tP,  the  exponent  now  positive.  Based  on  experience  at  Camp 
Century  and  the  Dye  Sites  in  Greenland,  I  expect  that  the  GLT  foundation 
will  experience  secondary  creep  as  the  main  mode  of  settlement  through¬ 
out  its  lifespan.  However,  during  dynamic  operations  (e.g.,  telescope 
pointing),  I  expect  that  the  snow  will  experience  some  amount  of  primary 
creep  as  it  receives  a  new  increment  of  load  during  rotation  of  the  eccentri¬ 
cally  loaded  dish.  As  such,  this  study  considers  models  that  capture  both 
primary  and  secondary  creep  behaviors. 


Figure  1.  Atypical  uniaxial  strain  curve  for  creep  loading  of  a  material. 


|< — - Secondary  - *\ - ► 

Primary  Tertiary 
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1.2  Objective 

The  objective  of  this  effort  was  to  develop  a  constitutive  model  that  de¬ 
scribes  the  creep  response  of  dense  snow  (p  =  440  -  890  kg/  m3)  for  use  in 
a  finite  element  analysis  to  quantify  long-term  settlement  of  a  structure 
placed  on  a  snow  foundation. 

1.3  Approach 

The  primary  creep  model  formulated  here  is  based  on  the  general  creep 
power  law.  To  model  secondary  creep,  this  effort  generalized  the  model  for 
snow  proposed  by  Mellor  and  Smith  ( 1966)  so  that  creep  could  be  com¬ 
puted  as  a  function  of  applied  load,  snow  temperature,  and  density  (i.e., 
following  the  form  of  Glen's  creep  law  for  ice,  Glen  1955) .  From  this,  I  for¬ 
mulated  a  general  model  that  is  the  basis  for  both  the  primary  and  second¬ 
ary  creep  models,  though  the  parameters  needed  for  each  differ.  Available 
laboratory  data  helped  to  determine  parameters  for  each  model.  The  form 
of  the  models  developed  in  this  study  are  enough  similar  to  the  generalized 
creep  model  used  in  the  ABAQUS  finite  element  software  that  the 
ABAQUS  creep  model  was  used  with  minimal  adaptation.  Finally,  I  ran 
simulations  of  test  cases  using  the  finite  element  method  to  validate  the 
models  against  available  data. 
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2  Primary  Creep 

Though  primary  creep  in  snow  has  not  been  extensively  studied,  a  few 
published  works  have  measured  and  developed  models  for  primary  creep 
of  snow.  Meussen  et  al.  ( 1999)  looked  at  the  creep  response  of  varying- 
density  snow  subjected  to  increasing  applied  loads  at  a  constant  tempera¬ 
ture  of  -20°C.  Chandel  et  al.  (2007)  provides  data  for  a  range  of  tempera¬ 
tures  (-3°C  to  -9°C)  while  holding  the  load  and  density  constant.  Combin¬ 
ing  these  data  sets  provides  an  opportunity  to  explore  the  variation  in 
strain  rate  during  initial  creep  loading  for  a  wide  range  of  density,  p ;  ap¬ 
plied  stress,  <t,  and  snow  temperature,  T.  Assuming  the  strain  rate  during 
primary  creep  will  have  the  following  dependency  e  -  /(t,  p,  a,  T),  I  pro¬ 
pose  the  use  of  the  power  form  creep  law  (ABAQUS  2003) : 

r  =  CpCT(jqtp  (1) 

I  chose  this  form  because  of  its  compatibility  with  the  standard  ABAQUS 
creep  model  (discussed  in  more  detail  in  the  section  "Implementation  in 
ABAQUS")  and  its  similarity  to  the  secondary  creep  model  proposed  by 
Mellor  and  Smith  ( 1966)  for  snow  and  Glen's  creep  law  for  ice  (more  com¬ 
monly  referred  to  as  Glen's  law) . 

Mellor  and  Smith  ( 1966)  proposed  that,  as  with  Glen's  law  for  ice,  the  Ar¬ 
rhenius  eguation  can  also  be  used  to  guantify  the  temperature  dependence 
of  creep  for  snow: 


k=  A  exp (-Q/RT) 


(2) 


where 


k  =  rate  constant  for  chemical  reactions, 

A  =  constant  or  proportionality, 

Q  =  activation  energy  (Cal/ mol), 

R  =  the  ideal  gas  constant  (1.985  Cal/  mol  K),  and 
T  =  absolute  temperature  in  Kelvin. 

Applying  this,  I  assume  that  the  temperature  coefficient  in  Eguation  ( 1) 
has  the  form  Ct  =  exp(-Q/  RT) .  Mellor  and  Smith  ( 1966)  also  proposed 
Q  =  f(/r,  cr)  where  p  is  the  snow  density  and  o  is  the  applied  stress.  Table  1 
presents  the  values  of  Q  determined  by  Mellor  and  Smith  ( 1996)  for  snow. 
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Here  I  make  the  assumption  that  Q  =  f (p)  only  and  that  the  variation  in  Q 
in  the  last  two  entries  in  Table  1  is  not  related  to  the  applied  stress  but  ra¬ 
ther  reflects  measurement  accuracy.  I  further  assume  that  the  average  Q 
for  p  =  832  kg/  m3  is  the  correct  value  for  that  snow  density.  This  adjust¬ 
ment  to  the  approach  of  Mellor  and  Smith  ( 1966)  makes  Ct  a  term  that  de¬ 
pends  on  snow  temperature  and  density  only.  With  this  adjustment,  I  per¬ 
formed  a  curve  fit  to  the  Q  data  in  Table  1  (the  details  of  which  are  de¬ 
scribed  in  the  next  section)  and  determined  empirically 

<2(kcal/mol)  =  1208.2 p-°-704;  p  (kg/m3).  (3) 

I  note,  that  applying  an  assumption  Q  =  constant  (e.g.,  average  of  all  Q  val¬ 
ues  in  Table  1)  degraded  the  model's  ability  to  replicate  measured  data, 
providing  a  strong  argument  that  Q  does  vary  with  snow  density  as  guanti- 
fied  in  Eguation  (3). 

Meussen  et  al.  ( 1999)  fit  their  strain  data  to  an  eguation  of  the  form 
s  =  Ctn.  The  strain  rate  is  guite  simply  the  time  derivative  of  this  expres¬ 
sion:  e  =  nCt11-1.  Therefore,  I  can  use  values  of  C  and  n  determined  by 
Meussen  et  al.  ( 1999)  to  get  an  estimate  for  Q,  and  p  =  n- 1  for  Eguation 
(1). 


Table  1.  Activation  energies  determined  by  Mellor  and  Smith  (1966). 


Snow  density,  kg/m3 

Applied  uniaxial  stress,  kPa 

Activation  energy,  kcal/mol 

436 

49 

17.8 

531 

49 

14.0 

644 

49 

13.4 

832 

98 

10.7 

832 

49 

11.9 

I  adjusted  the  fit  used  by  Meussen  et  al.  ( 1999)  to  more  of  a  creep  power 
law  form  (e  =  C'Ci  crri11)  and  recomputed  the  constant  after  accounting  for 
the  inclusion  of  stress  and  temperature.  Furthermore,  I  determined  that 
the  exponent  g  =  1  in  Eguation  ( 1)  describes  the  response  of  the  data  of 
Meussen  et  al.  ( 1999) .  Because  all  of  the  data  acguired  by  Meussen  et  al. 
(1999)  was  for  snow  at  -20°C,  I  included  data  from  Chandel  et  al.  (2007) 
for  snow  at  temperatures  of  -3°  C  to  -9°  C  and  obtained  values  of  C'  and  n 
for  these  data,  also.  This  analysis  reveals  that  C'  varies  with  the  snow  den¬ 
sity  according  to  the  empirical  relationship  C'  =  2.5  x  1014  exp(-0.046p); 
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R2  =  0 .98 .  Furthermore,  n  varies  only  slightly  with  snow  density  and  tem¬ 
perature;  the  average  value  is  n  =0.526  (standard  deviation  of  0.0834), 
and  I  assume  that  n  =  0.5.  The  strain  rate  was  determined  by  p  =  n  -  1  = 
-0.5  and  CP  =  nC'  =  1.25  x  10 14  exp(- 0.046/?),  and  thus  the  model  for  pri¬ 
mary  creep  of  snow  is 

£  —  CpCTotv\ p  -  —0.5;  Ct—  exp(—  Q/ RT)] 

CRs-^Pa-1)  =  1.25  x  1014  exp(-0.046/?)  (4) 

where  Q  is  determined  from  Eguation  (3) .  I  provide  a  comparison  of  the 
measured  strain  to  this  model  in  Figure  2.  The  data  are  nondimensional- 
ized  in  this  plot  to  =  t0  5  so  that  they  can  be  plotted  against  the  single 

line  (y  =  t0-5)  that  represents  the  model  response.  Note  that  there  is  con¬ 
siderable  scatter  in  the  data  and  that  the  collapse  is  not  as  tight  as  I  would 
hope  when  considering  the  fit  to  C';  and  by  extension  CP,  has  such  a  high 
coefficient  of  determination  (R2).  Yet,  clearly  the  exponent  p  =  -0.5  is  a 
reasonable  representation  of  the  trend  in  the  data.  The  scatter  in  the  data 
is  likely  a  result  of  differences  in  the  snow  structure  and  the  degree  of  in¬ 
terparticle  sintering,  which  is  the  type  of  variability  I  expect  in  the  natural 
environment;  and  I  consider  Eguation  (4)  adeguate  for  describing  general 
trends  of  the  primary  creep  response  of  snow.  However,  clearly  further 
work  is  necessary  to  account  for  the  sources  of  variability  evident  in  Figure 
2. 

Note  that  this  model  is  a  little  tenuous  as  the  data  I  was  able  to  locate  ex¬ 
tend  from  only  -3°C  to  -20°C.  I  would  like  to  be  able  to  extend  the  model 
to  give  us  an  estimate  of  the  primary  creep  behavior  around  -  30  °C  (the  ex¬ 
pected  temperature  at  depth  for  the  GLT  foundation) .  For  now,  I  restricted 
the  applicability  of  Eguation  (4)  to  -3°C  >  T  >  -50°C  with  the  understand¬ 
ing  that  this  model  reguires  updating  as  further  data  become  available. 
Even  considering  these  limits,  the  model  appears  to  perform  adeguately 
for  snow  densities  between  440  and  890  kg/  m3,  which  encompass  the 
range  reguired  for  this  study. 


Normalized  Axial  Strain 
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Figure  2.  Comparison  of  the  primary  creep  model  to  experimental  data  for  several  tempera¬ 
tures.  The  first  number  in  the  legend  is  temperature  (°C),  the  second  is  snow  density  (kg/m3), 
and  the  third  is  the  applied  pressure  (kPa).  (CSU  is  Chandel  et  al.  2007,  and  MMO  is  Meussen 

et  al.  1999). 
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3  Secondary  Creep 

Mellor  and  Smith  ( 1966)  proposed  the  Arrhenius  equation  to  quantify  the 
secondary  creep*  behavior  of  the  snow: 

s  =  Cxe  Q/rt  (5) 

Mellor  and  Smith  ( 1966)  empirically  determined  but  did  not  publish  the 
constant  Ci,  units  of  s~\  This  is  a  unique  application  of  the  Arrhenius 
equation;  rather  than  explicitly  account  for  the  effect  of  applied  stress  as  is 
done  in  the  creep  power  law  (Equation  [1]),  Mellor  and  Smith  ( 1966)  pro¬ 
posed  a  dependence  of  Q  on  applied  stress  as  discussed  above.  I  propose  a 
modification  to  the  approach  of  Mellor  and  Smith  ( 1966)  as  follows.  As 
with  the  primary  creep  model  I  explicitly  account  for  applied  stress  via 
Ci  =  Kcrto  yield  Glen's  law  with  the  exponent  on  stress,  q  =  1: 

£  —  Kae  Q/rt  (6) 

Further,  I  assume  that  K  is  a  function  of  only  snow  density.  Glen's  law  was 
developed  to  describe  the  creep  response  in  ice;  here,  I  assume  that  it  can 
also  be  applied  to  describe  the  secondary  creep  response  of  dense  snow. 

To  determine  how  the  parameter  K  varies  with  snow  density,  I  relied  on 
the  data  of  Mellor  and  Smith  ( 1966)  presented  in  Figure  3.  The  present 
study  determined  the  parameter  K  from  a  least  squares  fit  of  Equation  (6) 
by  using  the  values  of  Q  =f(/i)  and  o  from  Table  1  (noting  that  for  p  = 

832  kg/  m3,  Q  is  the  mean  value  for  the  two  stress  levels) . 

I  observed  that  the  data  for  K  and  Q  associated  with  a  snow  density  of 
531  kg/  m3  had  large  residuals  in  comparison  with  the  other  cases.  This  had 
the  effect  of  degrading  the  performance  of  the  model  when  this  study  used 
fits  obtained  for  K  and  Q  =  f \p)  such  that  the  model  did  not  faithfully  fol¬ 
low  any  of  the  cases  presented  in  Figure  3. 


*  Mellor  and  Smith  (1966)  note  that  it  took  strains  of  2xl0~2  to  10-1  to  transition  to  the  secondary  creep 
stage  in  their  uniaxial  compression  tests. 
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Figure  3.  The  symbols  represent  the  temperature  dependent  creep  data  of  Mellor  and  Smith 
(1966).  The  lines  indicate  the  predicted  value  from  Equation  (2).  The  applied  stress  for  all 

cases  is  49  kPa  except  where  noted. 


Temperature,  °C 


Therefore,  I  corrected  this  by  adjusting  the  values  of  K  and  Q  for  the  case  p 
=  531  kg/  m3  to  improve  the  fit  through  the  remaining  points.  This  greatly 
improved  the  model's  ability  to  reproduce  all  of  the  cases  given  in  Figure  3, 
except  for  p  =  531  kg/  m3.  The  resulting  model  is 


E  =  KCto ;  A"(s_1Pa_1)  =  3.32x10 77/?~2a3;  p  (kg/m3),  (7) 

which  is  very  similar  to  Eguation  ( 3),  with  p  =  0,  which  removes  the  time 
dependence  on  strain  rate  as  is  characteristic  of  secondary  creep.  Figure  4 
compares  the  model  results  to  the  data  for  snow  temperatures  of  -  1°C  to 
-35°C.  The  model  was  run  for  8  cases  {p  =  436,  500,  531,  600,  644,  700, 
800,  and  832  kg/  m3)  and  an  applied  stress  of  49  kPa.  Figure  4  shows  a 
smoothly  varying  model  result. 

In  developing  Eguation  (7),  I  assumed  that  the  model  should  be  a 
smoothly  varying  function  over  the  density  range.  This  was  in  contrast  to 
the  data  presented  in  Figure  3  that  show  a  large  drop  in  the  strain  rate  as¬ 
sociated  with  a  change  in  density  of  about  100  kg/  m3  between  p  =  436  and 
531  kg/  m3  and  a  very  small  drop  in  strain  rate  for  a  similar  magnitude 
density  change  from  531  to  644  kg/  m3.  Figure  4  also  shows  how  the  data 
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for  p  =  531  kg/  m3  departs  from  the  overall  trend  in  the  data.  This  indicates 
either  that  there  was  a  particular  problem  with  the  data  obtained  by  Mellor 
and  Smith  ( 1966)  associated  with  p  =  531  kg/  m3  (e.g.,  the  grain  structure 
of  the  synthesized  snow  particles  were  inconsistent  with  the  other  sam¬ 
ples)  or  that  there  is  an  abrupt  change  in  the  response  of  snow  in  the  range 
of  436  and  531  kg/  m3.  The  current  form  of  the  model  does  not  capture  that 
change  if  the  latter  is  true. 

Figure  4.  Comparison  of  data  of  secondary  creep  from  Mellor  and  Smith  (1966)  ( symbols )  to 
Equation  (7)  (line)  for  the  snow  temperatures  indicated  and  applied  stress  of  49  kPa.  Note 
that  for  the  temperature  of  -20°C,  the  ABAQUS  constitutive  model  that  follows  was  applied  to 
a  single  element  model  and  that  the  results  between  the  finite  element  model  and  Equation 
(7)  are  indistinguishable.  For  all  other  lines  in  this  figure  except  -20°C,  only  the  model  (Equa¬ 
tion  [7])  results  are  presented. 
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4  Implementation  in  ABAQUS 

ABAQUS  finite  element  software  provides  several  creep  models  for  simu¬ 
lating  the  time- dependent  response  of  materials.  The  time- dependent  law 
(ABAQUS  syntax:  *  creep  ,  law=t  ime  )  model  has  the  creep  power  law  form 
(ABAQUS  2003) 


£  —  Aaqtp  (8) 

where  t  is  time,  A  is  a  parameter  similar  to  C  used  in  Equation  (7),  and  q 
and  p  determine  the  power  law  response  for  stress  and  time.  Equation  (4) 
and  Equation  (7)  are  essentially  in  this  form  with  q  =  1  for  both  primary 
and  secondary  creep  and  p  =  -  0.5  for  primary  creep  while  p  =  0  for  sec¬ 
ondary  creep.  A  =  f(T ,p)  for  both,  though  the  factors  Q,  and  K  differ  for 
primary  and  secondary  creep. 

In  addition  to  the  creep  response  of  the  material,  the  elastic  properties  also 
need  to  be  specified.  Shapiro  et  al.  (1977)  provided  a  compilation  of  pub¬ 
lished  material  properties  of  snow.  They  showed  that  the  elastic  modulus, 
Y,  varies  mainly  with  snow  density,  though  there  is  considerable  scatter  in 
the  data.  One  would  expect  that  there  might  be  temperature  dependence, 
also,  though  that  is  not  readily  apparent;  if  such  dependence  exists,  it  is 
obscured  in  the  scatter  of  the  published  results  and  may  explain  the  fact 
that  temperature  was  not  plotted  as  an  additional  parameter.  From  these 
data,  I  extracted  the  trend  in  Y  vs.  p  and  provide  this  as  a  look-up  table  in 
ABAQUS  (Table  2). 

Table  2.  Dependency  of  the  elastic  modulus  on  snow  density  extracted 
from  data  compiled  in  Shapiro  et  al.  (1997).  The  data  extracted  follows 
more  closely  the  data  reported  by  Shapiro  et  al.  (1997)  that  was  ac¬ 
quired  from  low-strain-rate  measurements,  more  closely  akin  to  the 
foundation  loading  considered  for  the  application  of  this  model. 


Snow  density,  kg/m3 

Elastic  modulus,  MPa 

200 

1.379 

380 

98.20 

410 

146.9 

470 

303.3 

540 

633.6 

620 

1319.0 

770 

4165.0 

900 

9532.0 
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Shapiro  et  al.  (1997)  also  compiled  the  available  data  for  Poisson's  ratio,  v, 
as  a  function  of  snow  density.  The  range  in  v  for  snow  density  from  250  to 
850  kg/m3  is  0.2-  0.4  with  no  clear  trend  with  density.  Therefore,  I  se¬ 
lected  a  constant  value  of  r = 0.3  for  these  simulations,  which  is  roughly 
the  mean  in  the  data. 

4.1  Primary  creep 

This  is  determined  by  applying  Eguation  (8)  with  g  =  1,  p  =  -0.5,  and  A  = 
CpCT,  which  can  be  specified  as  a  function  of  temperature  for  a  specified 

snow  density.  Below  is  an  example  of  the  material  model  properties  for 
snow  with  a  density  of  550  kg/m3  and  a  temperature  of  -30°C.  This  has 
been  implemented  as  parametric  input  in  ABAQUS  so  that  the  values  of  A 
and  Young's  Modulus  are  computed  automatically  when  provided  the 
snow  density  and  temperature. 


* MATERIAL,  NAME=SnowPC 
* CREEP,  LAW=TIME 


A (m-s/kg) 

q 

P 

T  (deg  C 

1 . 43571390929e-ll, 

1., 

LO 

o 

1 

-50  . 

2 . 9042702382e-ll, 

1., 

LO 

o 

1 

LO 

1 

5 . 69999427094e-ll, 

1., 

LO 

o 

1 

-40  . 

1.08744605408e-10, 

1., 

LO 

o 

1 

-35. 

2.02020966228e-10, 

1., 

LO 

O 

1 

-30  . 

3.66048681342e-10, 

1., 

LO 

O 

1 

LO 

C4 

1 

6.47855264808e-10, 

1.  , 

LO 

O 

1 

-20  . 

1 . 12152018338e-09, 

1., 

LO 

O 

1 

LO 
\ — 1 

1 

1 . 90140333825e-09, 

1., 

LO 

o 

1 

o 
\ — 1 

1 

2 . 85960135075e-09, 

1., 

LO 

o 

1 

-6. 

*ELASTIC 
719275000.0,0.3 
^DENSITY 
550 . 0 

The  above  model  was  applied  in  a  single  element  model  with  a  static  load 
applied  to  the  top  of  the  element  as  a  fixed  mass.  Therefore,  gravity  (g  = 
-9.8066  m/  s2)  acting  on  the  mass  provided  the  axial  stress  on  the  sample 
and  replicated  static  creep  loading.  Initially,  the  load  was  applied  using  a 
*  static  step  with  an  amplitude  ramp  so  that  the  acceleration  of  gravity  in¬ 
creased  from  0  to  g  in  1  s.  This  allowed  the  sample  to  deflect  elastically. 
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The  load  was  held  constant  during  the  next  step  ( *  vi  sco)  to  allow  the  sam¬ 
ple  to  deform  viscoelastically  in  creep.  The  duration  of  this  step  was 
1000  s.  The  deformation  rate  during  this  second  load  step  was  used  to 
compute  the  strain  rate. 

This  study  ran  four  cases  to  demonstrate  the  primary  creep  model.  The 
first  was  a  comparison  to  a  measured  data  set  obtained  by  Meussen  et  al. 

( 1999)  and  presented  in  Figure  5;  the  agreement  with  this  test  case  is  very 
good.  The  next  section,  "Comparison  of  Primary  and  Secondary  Creep 
Models,"  discusses  the  other  three  cases. 

Figure  5.  Implemented  in  ABAQUS,  a  comparison  of  model  results  (Equation  [3])  to  the  creep 
behavior  of  snow  observed  by  Meussen  et  al.  (1999).  The  simulated  and  measured  data  are 
for  a  snow  temperature  of -20°C  and  a  density  of  450  kg/m3.  The  applied  load  is  220  kPa. 


1  10  100  1000  10000  100000  1000000 

Time  (s) 


4.2  Secondary  creep 

For  this  case,  I  used  in  Eguation  (8)  g  =  1,  p  =  0,  and  A  -  K CT .  In 
ABAQUS,  A  can  be  specified  as  a  function  of  temperature  for  a  specified 
snow  density. 

Using  the  above  information,  the  material  specification  in  ABAQUS  for 
snow  at  -20°C  and  a  density  of  436  kg/  m3  is  as  follows: 


ERDC/CRREL  TR-17-7 


13 


* MATERIAL,  NAME=SnowSC 
*CREEP ,  LAW=TIME 


A (m-s/kg) 

q 

P 

T  (deg 

2 . 47086400868e-14, 

1., 

0., 

-50  . 

5. 66466039457e-14, 

1., 

0., 

LO 

1 

1.25324020658e-13, 

1., 

0., 

-40  . 

2 . 68166686721e-13, 

1., 

0., 

-35. 

5.56134501594e-13, 

1., 

0., 

-30  . 

1 . 11990624221e-12, 

1., 

0., 

-25. 

2 . 19365137898e-12, 

1., 

0., 

-20  . 

4 . 18635707718e-12, 

1., 

0., 

-15. 

7 . 79530507393e-12, 

1., 

0., 

o 
\ — 1 

1 

1 . 41825869341e-ll, 

1., 

0., 

-5. 

2 . 25326892625e-ll, 

1., 

0., 

-1 . 

*ELASTIC 

214673333.333,0.3 

^DENSITY 

436.0 

This  was  applied  in  the  same  single  element  model  described  above  and 
was  used  for  validation  of  the  primary  creep  model,  which  replicated  the 
loads  imposed  by  Mellor  and  Smith  ( 1966) .  The  finite  element  model  re¬ 
sult  was  compared  against  the  data  from  Mellor  and  Smith  ( 1966)  at 
-20°C,  and  it  exactly  reproduced  the  model  (Eguation  [7])  shown  in  Fig¬ 
ure  4. 

4.3  Including  variation  of  density  with  depth 

To  include  the  variation  of  density  with  depth,  I  added  an  additional  field 
variable  to  the  material  models  shown  previously:  depth.  This  was  imple¬ 
mented  with  a  user-defined- field  subroutine  in  ABAQUS  (usdfld  .  f— Ap¬ 
pendix  A  provides  an  example) .  The  usdfld  .  f  subroutine  used  for  this  ef¬ 
fort  extracts  the  depth  in  the  model  and  applies  the  properties  in  a  tabular 
fashion  similar  to  what  is  shown  above  for  temperature.  The  variation  for 
density  and  depth  provided  for  in  the  table  is  as  follows. 

A  density  vs.  depth  relationship  is  reguired,  such  as  data  collected  by  Al¬ 
bert  (unpublished)  at  Summit,  Greenland.  Haehnel  and  Knuth  (2011)  de¬ 
termined  a  curve  fit  to  the  Albert  (unpublished)  data: 
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p  =  323  +  68.26d045;0<d<  120m  (9) 

p  g|)  =  912;  d>  120m 

For  the  ABAQUS  material  model,  I  formulate  a  linear  interpolation  of  den¬ 
sity  variation  with  depth  for  three  depths,  z  =  f( surface  density),  20  m,  and 
120  m.  Though  this  may  seem  coarse,  it  does  an  adequate  job  of  following 
the  depth-  density  trend  over  the  region  of  interest  ( ~3  m  to  50  m)  as 
demonstrated  in  Figure  6.  Note  that  the  starting  depth  is  not  z  =  0  m  be¬ 
cause  the  foundation  will  be  set  on  a  snow  surface  that  is  below  the  natural 
grade  (i.e.,  the  top  surface  will  be  excavated  down  to  some  depth,  e.g.,  3  m) 
and  the  foundation  will  be  set  on  that  surface.  Note  that  if  the  ABAQUS 
model  required  a  better  definition  in  the  density-  depth  relation,  one  could 
increase  the  number  of  depths  in  the  input  table  for  ABAQUS. 


Figure  6.  Comparison  of  depth  vs.  density  at  Summit  Station,  Greenland  (solid line),  and  lin¬ 
ear  approximation  using  3  points. 
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The  material  model  for  secondary  creep  that  varies  with  depth  is  then  as 
follows: 


**  Depth  varying  density 
* MATERIAL,  NAME=SnowSC_DV 
*CREEP ,  LAW=TIME,  DEPENDENCE  IE  S=1 


A (m-s/kg) 

n 

m 

T  (deg  C) 

Depth (m 

4 . 29290734498e-13, 

1., 

o.. 

-30 . 5, 

4 . 55145 

8 . 47036799953e-14, 

1., 

o.. 

-30 . 5, 

20 .01 

5 . 94998797247e-16, 

1., 

o.. 

LO 

O 

00 

1 

120 . 01 

*USER  DEFINED  FIELD 

*DEPVAR 

1, 

*ELASTIC,  DEPENDENCIES=1 

**  Youngs  Poisson  T (deg  C)  Depth (m) :  Summit,  Greenland 


67 . 54e6, 

0 . 3, 

LO 

O 

00 

1 

0  . 

98 . 20e6, 

0 . 3, 

-30 . 5, 

0 . 670 

14  6. 9e  6 , 

0 . 3, 

-30 . 5, 

1.714 

303 . 3e6, 

0 . 3, 

-30 . 5, 

5.500 

633 . 6e6, 

0 . 3, 

-30 . 5, 

13 . 07 

1319. e6. 

0 . 3, 

-30 . 5, 

26.25 

4165. e6. 

0 . 3, 

-30 . 5, 

65.11 

9532 ,e6. 

0.3, 

-30 . 5, 

114.8 

*DENSITY 
458 . 0 


The  material  model  for  primary  creep  is  not  shown  here  but  is  formulated 
similarly 
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5  Comparison  of  Primary  and  Secondary 
Creep  Models 

To  verify  that  the  primary  and  secondary  creep  models  behaved  as  ex¬ 
pected,  I  conducted  finite  element  simulations  of  axial  compression  on  a 
radially  unconfined  snow  sample  by  using  Equation  (4)  (primary  creep) 
and  Equation  (7)  (secondary  creep).  For  the  primary  creep  model,  I  also 
explored  the  response  of  the  snow  with  temperature.  Figure  7  provides  the 
results  of  these  simulations. 

Figure  7  shows  that  the  axial  strain  decreases  with  temperature  as  would 
be  expected  and  thus  demonstrates  that  the  primary  creep  model  pre¬ 
sented  in  Equation  (4)  behaves  appropriately.  Second,  I  compared  the  ax¬ 
ial  strain  that  the  secondary  creep  model  (Equation  [7])  would  predict  over 
the  same  elapsed  time.  For  this  comparison,  I  used  a  constant  snow  tem¬ 
perature  of  -30°C  (the  red  solid  and  dashed  lines  in  Figure  7). 

Figure  7.  Comparison  of  axial  strain  predicted  using  Equation  (3)  (primary  creep)  as  a  func¬ 
tion  of  snow  temperature.  For  comparison,  the  figure  also  shows  the  axial  strain  predicted  for 
the  secondary  creep  relationship  {dashed line)  (Equation  [6]).  The  snow  density  for  all  of 
these  simulations  is  550  kg/m3.  The  applied  stress  is  15  kPa. 
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Clearly,  the  axial  strain  during  secondary  creep  is  much  smaller  than  that 
experienced  during  primary  creep.  I  note  that  the  time  duration  of  these 
simulations  is  a  little  more  than  8  hours  (roughly  eguivalent  to  a  GLT 
pointing  operation) .  One  can  observe  that  the  slope  of  the  primary  creep 
curve  for  -30°C  (red  solid  line)  at  the  end  of  the  simulation  (around  8 
hours)  is  still  higher  than  the  slope  of  the  secondary  creep  line,  though  the 
primary  creep  model  is  clearly  trending  toward  the  shallower  slope  of  the 
secondary  creep  line.  At  about  8  hours  of  simulated  time,  the  strain  rate 
associated  with  primary  creep  is  about  20  times  larger  than  the  strain  rate 
that  the  snow  would  experience  during  the  secondary  creep  regime.  When 
I  applied  Eguation  (4)  to  longer  times,  I  found  that  for  a  snow  temperature 
of  -30°C,  it  would  take  about  500  hours  (21  days)  for  the  strain  rate  asso¬ 
ciated  with  primary  creep  to  decline  to  the  same  level  as  that  seen  for  sec¬ 
ondary  creep  (Eguation  [7])  of  snow  at  the  same  temperature.  Note  that 
the  time  to  transition  from  primary  to  secondary  creep  for  these  simula¬ 
tions  is  not  universal  but  is  highly  dependent  on  the  load  configuration. 

For  a  uniaxial  compression  test  geometry  (results  in  Figure  6),  the  model 
predicts  the  time  to  transition  from  primary  to  secondary  creep  response 
is  about  21  days.  For  a  radially  confined  case  (e.g.,  a  building  placed  on  a 
semi- infinite  snow  foundation),  the  time  to  transition  from  primary  to  sec¬ 
ondary  creep  is  likely  longer. 
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6  Model  Validation  Against  Field  Data 

The  model  validation  discussed  in  previous  sections  compares  the  model 
results  to  uniaxial  compression  tests  performed  in  the  laboratory.  I  now 
show  results  for  the  secondary  creep  model  compared  to  field  measure¬ 
ments  made  by  Reed  ( 1966)  wherein  a  wooden  raft  is  loaded  with  a  weight 
and  placed  on  the  snow  surface  of  the  Greenland  ice  sheet  at  Camp  Cen¬ 
tury.  The  settling  of  the  raft  was  monitored  over  the  course  of  about 
2  years.  Secondary  creep  was  the  dominant  deformation  mechanism  dur¬ 
ing  these  experiments.  Therefore,  I  use  the  depth- varying  secondary  creep 
model  (i.e.,  snowsc_Dv  in  Section  4.3)  in  these  simulations.  The  density  of 
the  snow  surface  on  which  the  rafts  were  placed  was  about  458  kg/  m3 
(Reed  1974)  on  a  surface  that  was  at  least  0 .46  m  below  the  original  snow 
surface.  Reed  ( 1966)  studied  several  rafts  of  varying  geometry  and  load¬ 
ings.  I  consider  only  the  tests  conducted  by  Reed  ( 1966)  that  used  square 
rafts  ranging  from  0.4572  to  1.372  m  on  a  side  and  loaded  at  47.8  kPa. 
Reed  (1966)  reported  that  the  average  snow  temperature  was  -20°C;  the 
temperature  at  the  3  m  depth  varied  from  - 17°  C  to  -24°C  throughout  the 
test  period  (approximately  700  days). 

At  Camp  Century,  the  density  vs.  depth  relationship  is  different  from  that 
at  Summit.  Benson  ( 1962)  provided  some  data  on  how  density  varies  with 
depth  in  the  region  near  Camp  Century.  I  extracted  the  data  from  Figure 
52  in  that  work.  A  reasonable  fit  to  the  data  is 

p  A  =  355.8  +  60.17d0,474;  0  <  d  <  100m.  (10) 

Note  that  this  is  only  applicable  down  to  a  100  m  depth.  Beyond  about 
100  m,  the  change  in  density  with  depth  progresses  more  slowly  as  it  ap¬ 
proaches  the  density  of  solid  ice.  Yet  Equation  (10)  provides  sufficient  in¬ 
formation  for  this  validation  effort  as  the  model  depth  is  only  to  about 
14  m. 

Also  of  note,  the  observations  of  Benson  (1962)  indicated  that  at  depths 
below  3.2  to  4  m,  the  temperature  is  fairly  uniform,  at  about  -25°C,  down 
to  at  least  8-  9  m;  so,  though  Reed  ( 1966)  observed  a  near- surface  average 
snow  temperature  of  -20°C,  the  temperature  in  the  column  of  snow  af¬ 
fected  by  the  surface  load  is  colder  and  can  range  down  to  -25°C.  For  this 
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reason,  I  ran  the  simulation  with  several  snow  temperatures  to  try  to  un¬ 
derstand  the  influence  of  snow  temperature  on  settlement  rate  and  how 
that  compares  to  Reed's  (1974)  observations. 

Figure  8  plots  the  settlement  rate  for  the  four  square  rafts  loaded  at 
47.8  kPa  reported  by  Reed  ( 1974)  (indicated  by  the  diamonds) .  The  solid 
line  indicates  a  least- squares  linear  fit  to  the  data. 


Figure  8.  Comparison  of  model  results  to  the  data  of  Reed  (1974).  This  is  the  same  data  ob¬ 
tained  by  Reed  (1966)  with  the  settlement  rates  adjusted  to  a  snow  density  of  458  kg/m3. 


Figure  9  shows  the  finite  element  model.  Because  of  symmetry,  I  used  a 
quarter  space  of  the  simulated  geometry.  The  snow  is  modeled  as  a  quarter 
cylinder,  and  one- quarter  of  the  square  raft  is  loading  the  model  at  the  cor¬ 
ner  of  the  sector.  The  stress  bulb  in  the  snow  is  dearly  visible  in  Figure  9. 
The  radius  of  the  domain  is  2.74  m.  The  outer  edge  of  the  domain  uses  in¬ 
finite  elements  to  project  the  domain  in  the  radial  direction.  The  height  of 
the  domain  is  13.7  m.  This  is  suffidently  tall  that  the  stress  bulb  does  not 
effectively  see  the  lower  boundary  of  the  domain. 

I  ran  simulations  at  three  snow  temperatures  (-20°C,  -25°C,  and  -30°C) 
and  for  raft  sizes  of  0.457,  0.914, 1.37,  and  1.83  m.  I  extended  the  raft  size 
beyond  that  used  by  Reed  1966  to  get  a  feel  for  how  the  settlement  rate 
trended  as  the  contact  area  continued  to  grow.  I  also  ran  an  additional  case 
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to  verify  the  computational  mesh  shown  in  Figure  9  was  sufficiently  fine. 
This  added  case  had  a  mesh  that  was  twice  as  fine  as  that  shown  in  Figure 
9;  the  snow  temperature  was  -30°C,  and  the  raft  size  was  1.37m.  Compar¬ 
ing  these  results  to  the  same  conditions  on  the  coarser  mesh,  I  found  no 
difference  in  the  computed  settlement  rate  and  therefore  concluded  that 
the  0 . 17  m  mesh  shown  in  Figure  9  was  sufficiently  fine  and  subsequently 
ran  all  of  the  cases  at  that  resolution. 


Figure  9.  The  model  domain  used  for  simulating  the  load  cases  outlined  in  Reed  (1966).  The 
raft  is  not  shown  to  allow  viewing  of  the  stress  state  of  the  snow  immediately  under  the  raft. 
The  load  case  is  quarter  space  for  the  1.37  m  raft  loaded  at  47.8  kPa. 


S.  Mises 
(Avg:  75%) 

+4.5006+04 
+4.1256+04 
+3.7506+04 
+33756+04 
+3.0006+04 
+2  ,G25e+04 
+22516+04 
+  1.9786+04 
+  1  5016+04 
+  1.1266+04 
+7.5106+03 
+3.7616+03 
+  1.1836+01 


I Kj 

m 


ODB:  Reed-54jlew.odb  AbaquslStandard  0.14-1  Thu  Apr  10  1429:50  EDT2015 


Step:  Step- 1 ,  Auy  raft  load  to  snow  (elastic  response) 

Increment  1:Step  Time 
Primary  Var  S.Mic-e-:- 

Deformed  Var:  H  Deformation  Scale  Factor:  +1  000e+00 


Figure  8  compares  the  simulated  results  to  Reed's  (1966, 1974)  data.  The 
agreement  between  the  model  with  the  snow  temperature  of  -25°C  and 
that  of  the  Reed  data  is  remarkably  good  when  considering  the  uncertainty 
and  spatial  and  temporal  variability  of  the  snow  temperature  in  the  field 
case.  The  model  also  faithfully  reproduces  the  field  observation  that  the  in¬ 
creasing  footing  area,  with  constant  pressure  at  the  base  of  the  footing,  in¬ 
creases  the  settlement  rate;  the  total  load  on  the  footing  is  increased,  lead¬ 
ing  to  a  higher  settlement  rate  even  though  the  pressure  under  the  footing 
is  constant. 
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7  Conclusions  and  Recommendations 

I  determined  parameters  and  empirical  expressions  to  allow  use  of  the 
ABAQUS  creep  constitutive  law  (a  creep  power  law.  Equation  [  8])  to 
model  the  primary  and  secondary  creep  response  of  a  dense  snow  {p  =  440 
-  890  kg/  m3) .  Owing  to  the  lack  of  data  that  spanned  from  the  primary  to 
secondary  creep  response  of  the  snow,  I  developed  two  separate  models 
from  the  available  data,  one  for  primary  creep  (Equation  [4])  and  another 
for  secondary  creep  (Equation  [7]). 

Prior  efforts  to  model  primary  and  secondary  creep  in  snow  resulted  in 
models  that  differ  from  the  creep  power  law  used  in  ABAQUS.  In  this  ef¬ 
fort,  I  reformulated  the  models  used  by  others  to  determine  creep  power 
law  parameters  that  could  be  used  in  the  standard  ABAQUS  creep  model. 

Comparing  model-predicted  strain  and  settlement  to  published  data 
shows  that  the  secondary  creep  model  effectively  reproduces  the  limited 
available  laboratory  and  field  observations  of  creep  deformation  of  dense 
snow.  This  study  could  perform  only  a  limited  validation  of  the  primary 
creep  model  because  of  limited  data  documenting  the  primary  creep  re¬ 
sponse  of  snow.  However,  Figure  7  compares  the  relative  magnitude  of 
strain  to  expect  during  primary  and  secondary  creep  phases.  These  results 
follow  expectations,  even  if  Equation  (4)  cannot  be  fully  validated  at  pre¬ 
sent. 

The  models  developed  here  (primary  creep.  Equation  [4],  and  secondary 
creep.  Equation  [  7] )  are  useful  in  illustrating  the  relative  deformation  that 
can  occur  during  the  two  creep  regimes.  However,  it  is  difficult  to  simulate 
long-term  behavior  that  captures  the  transition  from  primary  to  secondary 
creep  response  of  the  snow  by  having  them  as  two  separate  models.  I  rec¬ 
ommend  a  followup  effort  to  merge  these  two  models  into  a  single  (e.g.. 
Maxwell- Voigt)  model  that  describes  the  creep  response  of  snow  from  ini¬ 
tial  load  application  through  primary  creep  phase  and  the  transition  to 
secondary  creep  deformation.  Development  of  a  more  generalized  snow 
creep  model  will  require  considerably  more  effort  than  that  expended  to 
determine  the  parameters  for  the  piecewise  models  developed  for  this 
study. 
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A  separate  but  parallel  work  applied  the  models  developed  in  this  effort  to 
estimate  the  differential  foundation  settlement  of  the  proposed  GLT  dur¬ 
ing  an  8- hour  pointing  operation. 
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Appendix  A:  Density  vs.  Depth  Material  Model 

The  subroutine  usdfld  .  f  is  given  below.  Note  that  it  is  set  up  such  that  the 
vertical  direction  (elevation  or  depth)  is  associated  with  they- axis: 
coord  ( 2 )  as  is  used  in  the  ABAQUS  model  for  which  it  is  paired  with.  The 
variable  eiev  is  the  elevation  in  the  model  that  represents  the  top  of  the 
snow  surface.  The  below  code  specifies  elev  =  0,  that  is,  the  model  is  set 
up  such  that  the  top  of  the  snow  surface  is  at  Y  =  0 .  This  may  not  be  the 
case  in  all  models  and  will  need  to  be  adjusted  accordingly  in  this  subrou¬ 
tine  to  address  individual  simulations.  Also,  the  starting  density  of  the 
snow  (snow  density  at  the  surface)  needs  to  be  entered  in  for  the  variable 
rho.  If  the  model  is  placed  on  the  surface  of  the  snow  and  is  not  excavated 
down  to  a  particular  depth  before  the  foundation  is  placed,  then  this  is  the 
surface  density  of  the  snow.  Otherwise,  it  is  the  density  at  the  depth  of  the 
excavation.  In  the  case  below,  rho  =  4  5  8  (units  of  kg/  m3) .  The  subroutine 
then  computes  the  density- adjusted  depth,  ddepth,  by  using  the  appropri¬ 
ate  density  vs.  depth  relationship.  In  this  case,  I  use  the  density  vs.  depth 
relationship  for  Summit  Station,  Greenland,  given  in  Haehnel  and  Knuth 
(2011).  If  a  different  relationship  is  appropriate,  it  needs  to  be  entered  into 
the  subroutine  here,  replacing  the  existing  expression. 

USDFLD.f: 

SUBROUTINE  USDFLD (FIELD, STATEV, PNEWDT, DIRECT, T, CELENT, 

+TIME, DTIME, CMNAME, ORNAME, NFIELD, NSTATV, NOEL, NPT, LAYER, 

+KSPT, KSTEP , KINC, NDI , NSHR, COORD, JMAC, JMATYP , MATLAYO, 

+LACCFLA) 

C 

INCLUDE  ' ABA_PARAM . INC ' 

C 

CHARACTER*  8  0  CMNAME , ORNAME 

CHARACTER* 3  FLGRAY (15) 

DIMENSION  FIELD (NFIELD) , STATEV (NSTATV) , DIRECT (3, 3) , 

+T  (3,  3)  ,  TIME  (2) 

DIMENSION  ARRAY (15) , JARRAY (15) , JMAC (*) , JMATYP (*) , 

+COORD (*) 

real  z,  elev,  ddepth,  rho 

c  snow  density 
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rho  =  458 . 

c  elevation  in  model  of  top  of  snow  surface 
c  this  needs  to  be  set  according  to  the  model  geometry 
elev  =  0. 

c  density  depth  d  =  (  [rho-323] /68 . 26) **  (1/0 . 45)  for  Summit, 
Greenland 

ddepth  =  (  (rho-323. ) /68. 26) **  (1/0.45) 

c  model  set  up  such  that  Y  is  the  vertical  direction 
z=coord (2 ) 

c  field  to  pass  back  to  ABAQUS 

FIELD ( 1 ) =ddepth  +  elev  -  z 

RETURN 

END 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining 
the  data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for 
reducing  this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA 
22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display 
a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


2.  REPORT  TYPE 


1.  REPORT  DATE  (DD-MM-YYYY) 
April  2017 


4.  TITLE  AND  SUBTITLE 


Technical  Report/Final 


A  Creep  Model  for  High-Density  Snow 


3.  DATES  COVERED  (From  -  To) 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 

Robert  B.  Haehnel 


5d.  PROJECT  NUMBER 

SAO-GLT-16-Ol 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Engineer  Research  and  Development  Center  (ERDC) 
Cold  Regions  Research  and  Engineering  Laboratory  (CRREL) 
72  Lyme  Road 
Hanover,  NH  03755-1290 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 

ERDC/CRREL  TR-17-7 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Smithsonian  Astrophysical  Observatory 
60  Garden  Street 
Cambridge,  MA  02138 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 

SAO 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 


14.  ABSTRACT 

This  effort  modifies  Mellor  and  Smith’s  creep  model  for  dense  snow  to  conform  to  the  more  general  creep  power  law  form 
(Glen’s  creep  law  for  ice  is  a  special  case  of  this  power  law).  The  present  study  used  this  general  form  as  the  basis  for 
developing  two  creep  models:  one  to  describe  the  primary  creep  and  another  that  models  secondary  creep  response  as  a 
function  of  applied  stress,  snow  density,  and  temperature.  Available  laboratory  data  helped  to  determine  the  parameters  for 
these  models.  These  models  were  recast  into  a  form  compatible  with  the  ABAQUS  finite  element  software  to  allow 
simulation  of  settlement  of  a  snow  foundation  subjected  to  a  long-term  load. 

The  secondary  creep  model  was  validated  against  the  observed  settlement  of  a  spread  foundation  on  snow  at  Camp  Century, 
Greenland.  Though  there  were  not  suitable  field  data  to  validate  the  primary  creep  model,  comparison  of  that  model  to 
results  obtained  for  the  secondary  creep  model  demonstrate  that  the  primary  creep  model  behaved  as  expected.  Additional 
work  is  necessary  to  merge  the  primary  and  secondary  creep  models  and  to  further  validate  the  primary  creep  response  of 
the  model  presented  here. 


15.  SUBJECT  TERMS 

Creep  deformation.  Finite  element  method.  Foundations— Cold  weather  conditions,  Greenland,  Polar  construction. 
Snow,  Snow  mechanics 


16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 

Unclassified 


b.  ABSTRACT 

Unclassified 


c.  THIS  PAGE 

Unclassified 


17.  LIMITATION 

18.  NUMBER 

OF  ABSTRACT 

OF  PAGES 

SAR 

34 

19b.  TELEPHONE  NUMBER  (include 
area  code) 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  239.18 


